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Abstract 



We provide a priori error estimates for variational approximations of the 
ground state energy, eigenvalue and eigenvector of nonlinear elliptic eigenvalue 
problems of the form — div(AVu) + Vu + f{u^)u = Xu, ||ii||i2 — 1. We focus 
in particular on the Fourier spectral approximation (for periodic problems) 
and on the Pi and P2 finite-element discretizations. Denoting by {us,Xs) a 
variational approximation of the ground state eigenpair (m, A), we are inter- 
ested in the convergence rates of \\us — ^H^fi, \\us — u||£_2, \Xs — X\, and the 
ground state energy, when the discretization parameter 5 goes to zero. We 
prove in particular that if A, V and / satisfy certain conditions, \Xs — X\ goes 
to zero as \\us — + Wu-s — ^(11^2. We also show that under more restric- 

tive assumptions on A, V and /, \Xs — X\ converges to zero as \\us — w||^i, 
thus recovering a standard result for linear elliptic eigenvalue problems. For 
the latter analysis, we make use of estimates of the error us — u in negative 
Sobolev norms. 



1 Introduction 

Many mathematical models in science and engineering give rise to nonlinear eigen- 
value problems. Let us mention for instance the calculation of the vibration modes of 
a mechanical structure in the framework of nonlinear elasticity, the Gross-Pitacvskii 
equation describing the steady states of Bose-Einstein condensates [9] , or the Hartree- 
Fock and Kohn-Sham equations used to calculate ground state electronic structures 
of molecular systems in quantum chemistry and materials science (see [3] for a 
mathematical introduction). 

The numerical analysis of linear eigenvalue problems has been thoroughly studied 
in the past decades (see e.g. [1]). On the other hand, only a few results on nonlinear 
eigenvalue problems have been published so far [13, 14]. 

In this article, we focus on a particular class of nonlinear eigenvalue problems arising 
in the study of variational models of the form 
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where 

is a regular bounded domain or a rectangular brick of M'^ and X = ifg (f^) 

or 

il is the unit cell of a periodic lattice TZ of M** and X = H^{n) 
with d = 1, 2 or 3, and where the energy functional E is of the form 

Eiv) = ^a{v,v) + ^ I F{v\x))dx 



n 



with 



a{u,v) = / {AVu) ■Vv+ Vuv. 
Jn Jn 



m Jn 

Recall that if is the unit cell of a periodic lattice TZ of M'', then for all s G M and 

fc G N, 

H^{n) = {v\n, V G H^l^iR'') I V 7^-periodic} , 
C^{n) = {v\n, V G C'^iK.'^) I V 7^-pcriodic} . 

We assume in addition that 

• A G (L°°(fi))''^'' and A(a;) is symmetric for almost all x G 17 (2) 
3a > s.t. C'^A(x)^ > a|^p for aU ^ G K'* and almost aU a; G O (3) 

• V G LP{n) for some p > max(l, d/2) (4) 

• G Ci([0,+oo),M) nC2((0,oo),R) and F" > on (0,+oo) (5) 
30 < g < 2, 3C G R+ s.t. Vi > 0, \F'{t)\ < C(l + 1«) (6) 
F"{t)t remains bounded in the vicinity of 0. (7) 

To establish some of our results, we will also need to make the additional assumption 
that there exists 1 < r < 2 and < s < 5 — ?' such that 

Vi? > 0, 3Ch G R+ s.t. VO < < i?, Vt2 e M, 

\F'{tl)t2 - F'{tl)t2 - 2F"{tl)t\{h -h)\< Cr (1 + \h - iiT- (8) 

Note that for aU 1 < m < 3 and all c > 0, the function F{t) = ci™ satisfies (5)- (7) 
and (8), for some 1 < r < 2. It satisfies (8) with ?- = 2 if 3/2 < m < 3. This allows 
us to handle the Thomas- Fermi kinetic energy functional {m = |) as well as the 
repulsive interaction in Bose-Einstein condensates (m = 2). 

Remark 1 Assumption ( 6) is sharp for d = 3, but is useless for d = 1 and can be 
replaced with the weaker assumption that there exist q < oo and C G M+ such that 
\F'{t)\ < C(l + t«) for all t G M+, for d = 2. Likewise, the condition < s < 5 - r 
in assumption (8) is sharp for d = 3 but can be replaced with Q < s < oo if d = \ 
ord=2. 

In order to simplify the notation, we denote by f{t) = F'{t). 

Making the change of variable p ^ and noticing that a(|ij|, \v\) = a{v,v) for all 
w G X, it is easy to check that 



/ = inf|£(p), p>0, ^p<^X, L^^^] 
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where 



We win see that under assumptions (2)-(6); (9) has a unique solution pQ and (1) 
has exactly two solutions: u = y/po and —u. Moreover, E is on X and for all 

V € X, E'{v) = AyV where 

A„ =-div(^V-) + l^ + /(w'). 

Note that defines a self-adjoint operator on L^(r2), with form domain X. The 
function u therefore is solution to the Euler equation 

Vu e X, {Auu - Xu, v)x',x = (10) 

for some A e R (the Lagrange multiplier of the constraint ||u|j|2 = 1) and equa- 
tion (10), complemented with the constraint ||u||l2 = 1, takes the form of the 
nonlinear eigenvalue problem 

AuU = \u , . 

In addition, u e C°(ri), u > Q in fl and A is the ground state eigenvalue of the 
linear operator A-^. An important result is that A is a simple eigenvalue of Au- It 
is interesting to note that A is also the ground state eigenvalue of the nonlinear 
eigenvalue problem 

search (/i, w) e M x X such that 

A.i,v = fiv (12) 
II^'IIl^ = 1, 

in the following sense: if {p., v) is solution to (12) then either p > X ot p, = X and 

V = ztu. All these properties, except maybe the last one, are classical. For the sake 
of completeness, their proofs are however given in the Appendix. 

Let us now turn to the main topic of this article, namely the derivation of a priori 
error estimates for variational approximations of the ground state eigenpair (A, u). 
We denote by (Xs)s>q a family of finite-dimensional subspaces of X such that 

min {\\u - vs\\h^, vs & Xs} — > (13) 

S i'O"^ 



(14) 



and consider the variational approximation of (1) consisting in solving 
Is = inf I^E{vs), vs ^ ^s, J ^5 = l| • 

Problem (14) has at least one minimizer ug, which satisfies 

Vi's e Xs, (Au,M5 - Xsus,vs)x',x = (15) 

for some Xg G M. Obviously, —us also is a minimizer associated with the same 
eigenvalue Xs- On the other hand, it is not known whether us and —us are the 
only minimizers of (14). One of the reasons why the argument used in the infinite- 
dimensional setting cannot be transposed to the discrete case is that the set 

{p I 3us e Xs s.t. \\us\\l^ = 1, p = ul} 

is not convex in general. We will see however (cf. Theorem 1) that for any family 
(mi5)i5>o of global minimizers of (14) such that {u, us) > for all S > 0, the following 
holds true 

\\us - uWffi — > 0. 

(5— j'0+ 
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In addition, a simple calculation leads to 

Xs - {{An - X){us - u),{us - u))x',X + 'Wu,us{u6-u) (16) 

Ja 

where 

Wu.us = Ug . 

US ~ u 

The first term of the right-hand side of (16) is nonnegative and goes to zero as 
\\us — It j 1^1. We will prove in Theorem 1 that the second term goes to zero at 
least as \\us — u||i6/(5-2,) . Therefore, jA^ — A| converges to zero with 5 at least as 

ll".? - + \\US - M||i6/(5-2,). 

The purpose of this article is to provide more precise a priori error bounds on 
\\s — A|, as well as on \\us — u\\fji, \\us — u\\i^2 and E(us) — E{u). In Section 2, 
we prove a series of estimates valid in the general framework described above. We 
then turn to more specific examples, where the analysis can be pushed further. In 
Section 3, we concentrate on the discretization of problem (1) with 



O = (0,27r)^ 

E{v) ^ \ [ |V«|2 + \ [ Vv^ + \ [ F{v^), 



in Fourier modes. In Section 4, we deal with the Pi and P2 finite element discretiza- 
tions of problem (1) with 

ri rectangular brick of M'', 
Lastly, we discuss the issue of numerical integration in Section 5. 



2 Basic error analysis 

The aim of this section is to establish error bounds on \\us — \\us — 

\Xs — X\ and E{us) — E{u), in a general framework. In the whole section, we make 
the assumptions (2)-(7) and (13), and we denote by u the unique positive solution 
of (1) and by us a minimizcr of the discretized problem (14) such that (us, u)]^2 > 0. 
We also introduce the bilinear form E"{u) defined on X x X by 

{E"{u)v,w)x',x = {AuV,w)x'.x + 2 / f'{u^)u^vw. 

When F G C2([0, -l-oo), M), then E is twice differentiable on X and E"{u) is the 
second derivative of £' at w. 

Lemma 1 There exists (3 > and M e such that for all v G X , 

< {{Au-X)v,v)x',x<M\\v\\]J^ (17) 
(i\\v\\]ii < {{E"{u)-X)v,v)x'^x <M\\v\\j,,. (18) 

There exists 7 > such that for all S > 0, 

7||it5 - u\\%i < {{Au - X){us - u), {us - u))x',x- (19) 
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Proof Wc have for all u G X, 



{{Au - X)v,v)x'.x < \\A\\l^\\Vv\\1, + \\V\\Ml..' + \\fiu')\\L^\\v\\ 



where p' ~ (1 — p ^) ^ and 

{{E"{u) - X)v,v)x'.x < {{Au - X)v,v)x'.x + 2||/'(ii2)^2||^^ ||^,||2 

Hence the upper bounds in (17) and (18). We now use the fact that A, the lowest 
eigenvalue of Au, is simple (see Lemma 2 in the Appendix). This implies that there 
exists 77 > such that 

yveX, {{Au-X)v,v)x',x>vi\\v\\l2-\{u,v)L2\^)>0. (20) 

This provides on the one hand the lower bound (17), and leads on the other hand 
to the inequality 



V« e X, {{E"{u) - X)v, v)x',x > 2 / f'{u^)u^v 



^ '/„,2^„,2„,2 

n 

As /' = F" > in (0, +00) and u > in il, we therefore have 
yveX\{0}, {{E"{u)-X)v,v)x',x>0. 

Reasoning by contradiction, we deduce from the above inequality and the first 
inequality in (20) that there exists ^ > such that 

yv e X, {{E"{u) - X)v, v)x',x > (21) 

Besides, there exists a constant C E M.+ such that 

yveX, {{Au-X)v,v)x'.x>^\\^v\\l,-C\\v\\l.. (22) 

Let us establish this inequality for d — 3 (the case when d = 1 is straightforward 
and the case when d = 2 can be dealt with in the same way). For all x Cz X, 

((A, ~ X)v, v)x',x = I {AVv) ■Vv+ f {V + f{v^) - X^ 

> c.\\Vv\\h-\\V\\L.\\v\\l.,,+if{0)-X)\\v\\l. 

> a\\Vv\\l. - \\Vh4v\\l~.'^lvf/r + (/(O) - X)\\v\\l. 

> a\\Vv\\h Cl/'\\V\\Ml'.'^lv\\'^r + (/(O) - X)\\v\\l, 



3-2p fscimff 



3/(2p-3) 



a 



where Ce is the Sobolev constant such that Vu £ X, \\v\\i^6 < CgWvWh'-- The coer- 
civity of E" {u) — A (i.e. the lower bound in (18)) is a straightforward consequence 
of (21) and (22). 

To prove (19), we notice that 

ll^alli^ - \{u,us)lA'^ > 1 - {u,us)l^ = ^\\us - u\\l2. 
It therefore readily follows from (20) that 

{{Au ~ X){us - u), {us - u))x'.x > - u\\l2. 

Combining with (22), we finally obtain (19). □ 
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For w G X' , we denote by ipw the unique solution to the adjoint problem 

find G such that , , 

yveu^, {{E"iu)-X)Ao,v)x',x = {w,v)x'^x, ^ ' 

where 

u-L = ji; G X I J uv = 

The existence and uniqueness of the solution to (23) is a straightforward consequence 
of (18) and the Lax-Milgram lemma. Besides, 

yw G L^n), < l3-'M\\w\\x' < I3-HI\\w\\l2. (24) 

We can now state the main result of this section. 



Theorem 1 Under assumptions (2)-(6) and (13), it holds 

\\us - uWhi — > 0. 

(5^0+ 

// in addition, (7) is satisfied, then there exists C G IR-|_ such that for all S > 0, 

'^\\us - u\\jji < E{us) - E{u) < ^\\us - u\\jji + C\\us - u||i6/(5-2,) , (25) 

and 

\Xs - A| < C {\\us - u\\l, + \\us - u|li6/(5-2,)) . (26) 

Besides, if assumption (8) is satisfied for some 1 < r < 2 and < s < 5 — r , then 
there exists (5o > and C G M+ such that for all < S < Sq , 

\\us-u\\hi < C min \\vs - u\\h^ (27) 



\u&-u\li2 < C (^\\us - u\\l2\\us - u\\ler/(5-s) 



+ 11^5 - fi||_H-i min \\iPus-u-i1js\\h^ ]■ (28) 

Lastly, if F" is bounded in the vicinity of 0, there exists C G M+ such that for all 
6>0, 

^Wus - u\\l, < E{us) - E(u) < C\\us - u\\l,. (29) 



Remark 2 //0<r + s<3, then 

\\u5 - w||26./(5-=) < \\us - u\\^l^^''^''^^^\\us - u\\^^e^^^'^^^ < \\us - u\\l2\\us - 
so that (28) implies the simpler inequality 

\\us - u\\\2 < C\\us - u\\hi min Wipus-u - i'sWn^ ■ (30) 
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Proof of Theorem 1 We have 

E{us)-E{u) = '^{AuUs,us)x',x -]^{AuU,u)x'^x 

F [us) - F {u') ~ / {u') {ul - u^). (31) 



n 



2 

Using (19) and the convexity of F, we get 

E{us)-E{u)>^-\\us-u\\]j,. 

Let lisu E Xs be such that 

\\u - UsuWhi = min{||w - vs\\hi,vs G Xs} . 

We deduce from (13) that (nsu)s>o converges to u in X when 6 goes to zero. 
Denoting by us = \\Ilsu\\J^llIsu (which is weh defined, at least for S small enough), 
we also have 

lim llu^ — ullffi = 0. 

(5^0+ 

The functional E being strongly continuous on X, we obtain 

\\us - u\\l^ < - {E{us) - E{u)) < - {E{us) - E{u)) 0. 
It follows that there exists Si > such that 

VO < 5 < (5i, \\us\\m <2\\u\\hi, \\us-u\\hi<^- 

We then easily deduce from (31) the upper bounds in (25) and (29). 
Next, we remark that 

Aa - A = (E'{us),us)x'.x ~ {E'{u),u)x',x 

= a{us,us) ~ a{u,u) + / f{ul)u]- / f{u^)u^ 

Jn Jn 

= a{us - u,us ~ u) + 2a{u,us - u) + / J{uf}u'l- / f{v?)v? 

Jn Jn 

= a{us — u,us — u) + 2X / u{us — u) — 2 / f{u^)u{us — u) 

Jn Jn 



Jn Jn 

= a{us — u,us ~ u) — X\\us — u\\ j^2 — 2 j f {u'^)u{us — u) 



n 



+ / f{ui)ui / f{u')u' 
Jn Jn 

= {{Au - \){us - u), {us - u))x',x + I Wu.us {us - u) (32) 



n 



where 



, f{uj)~f{u') 

Wu.us = Us . 

Us — u 



As u G L°°(il), wc have 



^u,us I — 



12u sup F"{t)t if \us\ < 2u 

t6(0:4||«||2^] 

2(1/(^2)1+ max 1/(01 I kil ii\us\>2u, 



and wc deduce from assumptions (6)- (7) that 

\Wu,uA<C{l+\us\'"'+'), 

for some constant C independent of 6. Using (17), we therefore obtain that for all 
< (5 < 

IA5-AI < A/|1U5 - + |lm„.„J|^6/(2,+ l)||M5 - u||^6/(5-2,) 

< M\\us - u\\]ji + C(l + \\us\\^^t^)\\us - u||2,6/(5-2,) 

< C - ujl^i + - ui|i,6/(5-2,)) , (33) 

where C denotes constants independent of 5. 

In order to evaluate the ^-norm of the error ug — u, we first notice that 

Vui e \\us~u\\h^ <\\us-vs\\h^ + \\vs~u\\h^, (34) 

and that 

\\us-vs\\]ji < I3~^ ({E"{u) - X){us - vs),{us - vs))x',x 
= [{{E"{u) - X){us - u), [us - vs))x',x 

+ {{E"{:u) - X){u - vs), {us - vs))x'^x] ■ (35) 



For all ws E Xs 
{{E"{u) - \){us - u),ws)x',x 

= - / {f{us)us - ,f{u'^)us - 2f'{u'^)u'^{u5 - u)) ws + {\s - X) / usws- (36) 
Jn Jn 

On the other hand, we have for all vs € Xs such that ||iJ5j|^2 = 1, 
us{us - vs) ^1 - I usi's = l-\\us - vs\\'i2. 



n Jn 



2' 



Using (8) and (33), we therefore obtain that for all < 5 < Si and all vs G Xs such 
that \\vs\\l^ = 1, 



\{{E"{u) - X){us - u), {us - vs))x',x\ < c(^\\us - u||26,v(5-=) - vs\\ 
+ {\\us - u\\jfi + \\us - u||i6/(5-25)) \\us - vs\\j^2^. 



(37) 



It then follows from (18), (35) and (37) that for all < 5 < 61 and all vs € Xs such 
that \\vs\\l^- = 1, 

\\us - vsWh^ < C {\\us - ull^i + \\us - u\\hAW& - ""sWm + hs - u\\m) ■ 
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Combining with (34) wc obtain that there exists Q < 82 < Si and C E M+ such that 
for all < (5 < (52 and all vs G Xg such that = 1, 

\\us - -ullffi < C\\vs - uWh^- 

Hence, for all < (5 < ^2 

\\us — u\\ffi < CJs where Js = min \\vs — wH/fi- 

vs<^Xs I ||uj|| 2 = 1 



We now denote by 



Js = min \\vs - u\\m, 



and by a minimizer of the above minimization problem. We know from (13) that 
Ug converges to u in when 5 goes to zero. Besides, 

Js < H/\\u''s\\l2 -u\\hi 

< Ilu?-1tllt7 



< 



< 



-'-111 nil 







H\ 


L2 


H\ 




H\ 




1 Js- 





For < S < 62 < 5i, we have — < \\us — < 1/2, and therefore 

< +1/2 and |lit!^||L2 > 1/2, yielding Js < 2{\\u\\hi +l)Js- Thus (27) 

is proved. 

Let Ug be the orthogonal projection, for the inner product, of us on the affine 
space {v e L^{Q) \ J^uv ^ 1^. One has 

u*seX, ug-ueu-^, u*s - us ^ -^\\us - u\\j^2U, 

from which wc infer that 

Wus-uWlo. = I {us - u){u*s - u) + I {us - u){us - Us) 
n 

{us - u){u*s - w) - -\\us - u||^2 / {us - u)u 
{us - u){u*s - u) + - ?i|||2 ^ ^ 

{us - u){u*s -u) + ^\\us - u\\l2 
= {us - u,u*g - u)x',x + ^\\us - u\\\2 
= {{E"{u) - X)ipus-u,u*s - u)x',x + ^\\us - u\\l2 



= {{E"{u) - X){us - u),1pus-u)x',X 

+ - u\\l2{{E"{u) - \)u,'iIju6-u)x\x + ^\\us - u||i2 
= {{E"{u) - X){us - u),ilju/i-u)x',x 



1, 



+ \\us-u\\^2 J f {u )u + - u|1l2. 
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For all -05 G Xs, it therefore holds 

\\us-u\\'l2 = {{E"{u) - \){us - u),ijjs)x',x 

+ {{E"{u) - X){us - u),ijjus~u - '^s)x',x 

+ \\us - u\\\2 I f{u^)u^ipus-u + ^\\us - 



n 4' 



L2- 



From (36); we obtain that for all tps G Xs n u^, 

{{E"{u)~ X){us-u),yjs)x'^x 

= - f {f{uj)us - f(u^)us - 2f'{u^)u^{us - u)) + (A^- - A) / {us - u)^s 
Jn Jn 

and therefore that for all t/js E Xs O , 

\{{E"(u)-X){us~u),i;s)x',x\ < C (^\\us - u\\l,^,,,_^, 

+ \\us - u||l6/5 {\\us - u\\jji + \\us - ui|2^6/(5-2,)) ^ llV'alUi- (38) 

Let ipg € XsH u-^ be such that 

Hus-u- iPsWh^ ^ min I] - "05 ll^i . 

4>s£Xsnu^ 

Noticing that ||V'^||hi < || V'ua-ullHi — P~^M\\us — we obtain from (18) and 

(38) that there exists C G M_|- such that for all < (5 < (Si, 



\\us-u\\l2 < C'\^\\us - u\\l2 

X [\\US - u\\ler-/(5-s) + \\US - u||l6/5 {\\us - + \\US - 'u||i6/(5-2,))) 

+ \\us - u\\hA\'4'us-u - V'°I|hi + \\us - uWh + \\us - u\\'^2 
Therefore, there exists < Sq < 62 and C G such that for all < S < Sq, 

\\US - u\\l2 < C (^\\us - u\\l2\\us - u\\ler/^5-s) + \\us - u\\hi\H'us-u - 

Lastly, denoting by ^Xs ^^'^ orthogonal projector on Xs for the inner product, 
a simple calculation leads to 

, / |in5fM|lHi\ 

Vu G u , mm \\vs - vWhi < U + htto" 112 ™5 \\v5-v\\m, (39) 

vs^Xsnu^ \ \\lLj^^U\\j^2 I veeXs 

which completes the proof of Theorem 1. □ 

Remark 3 In the proof of Theorem 1, we have obtained bounds on \\s — A| from 
(32), using LP estimates on Wu^us o.'^d {ug — u) to control the second term of the 
right hand side. Remarking that 



f{u^)u - f{u l)u^2f'{ ul)ul{u - us) 
{us - uY 



^Wu,us = -U^^ \; -^US 



f{ul)us - f{u^)us - 2f'{u^)u^us -u) ^ 

-Us -. Vu 

{us - uy 

+2uus (f'{u'j) Vus + f'{u^)Vu) + 2us ~ -^^^ ^ Vu^ . 

Us — u 
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we can see that if us is uniformly bounded in L°°{Cl) and if F satisfies (8) for r = 2 
and is such that F"(t)t^/'^ is bounded in the vicinity ofO, then Wu.us is uniformly 
bounded in X . It then follows from (32) that 

\\s ~\\<C {\\us - + \\us - u\\x') , 

an estimate which is an improvement of (26). In the next two sections, we will see 
that this approach ( or analogous strategies making use of negative Sobolev norms of 
higher orders), can be used in certain cases to obtain optimal estimates on \Xs — A| 
of the form 

\\5-X\<C\\us-u\\]j^, 
similar to what is obtained for the linear eigenvalue problem — Au + Vu ~ Xu. 



3 Fourier expansion 

In this section, wc consider the problem 

M I^E{v), V e X, J^v^ ^1^ , (40) 

where 

= (0, 27r)'^, with d= 1, 2 or 3, 

X = H^^{n), 

We assume that V G II^{il) for some cr > d/2 and that the function F satisfies (5)- 
(7), (8) for some 1 < r < 2 and < s < 5 - r, and is in cM+i^'"~M+'((0, +oo),M) 
(with the convention that C'''" = C'' if A: e N). 

The positive solution u to (40), which satisfies the elliptic equation 

-Au + Vu + f{u^)u = Am, 

then is in H^'^^ifl) and is bounded away from 0. To obtain this result, we have 
used the fact [11] that if s > d/2, g e CW'''-[''1+"(R, R) and v e H^{n), then 
g{v) e H^ifl). 

A natural discretization of (40) consists in using a Fourier basis. Denoting by 
ek{x) = (27r)-^/2e''=-^, wc have for aU v e L^ifl), 

vix) = ^ Vkek{x), 



where Vk is the k Fourier coefficient of v: 

Vk:= / v{x) ek{x)dx = {2tt)-^/'^ [ v{x)i 

The approximation of the solution to (40) by the spectral Fourier approximation is 
based on the choice 

Xs^Xn ^ Span{efc, < N}, 
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where |fc|* denotes either the /^-norm or the /°°-norm of k (i.e. either |fc| = 
(^i^i\ki\'^y^^ or |fc|oo = niaxi<i<(i |fci|). For convenience, the discretization pa- 
rameter for this approximation will be denoted as N. 

Endowing _ff^(ri) with the norm defined by 

1/2 

hh^ = [ e {^+\k\iy\vy 

we obtain that for all s g M, and all v G H^{n), the best approximation of v in 
H^{^1) for any r < s is 

keZ'',\k\,<N 

The more regular v (the regularity being measured in terms of the Sobolev norms 
If), the faster the convergence of this truncated series to v: for all real numbers r 
and s with r < s, we have 

V«ei/^(r!), \\v-Il^v\\Hr<j;^\\v\\H^. (41) 



Let ujv be a solution to the variational problem 
inf <^ E{vn), vn G Xn, 



— 

in 

such that (ujv,u)i2 > 0. Using (41). we obtain 

\\u-UnU\\h1 < ^l^-^ \\u\\H<r + 2, 

and it therefore follows from the first assertion of Theorem 1 that 

lim \\un — u\\hi = 0. 

N^oc 

We then observe that ujv is solution to the elliptic equation 

-Aupf + Un [Vun + f{u%)uN] = XnUn- (42) 
Thus Mat is uniformly bounded in iJ^(f7), hence in L°°{il), and 

A{un-u) = Un {V{uN - u) + f{u%)uN - f{u'^)u) 

-(/ - UnKVu + f{u^)u) - \n{un -u)- {Xn - X)u. (43) 

As {uN)NeTi in bounded in L°°(f2) and converges to u in H^{fl), the right hand 
side of the above equality converges to in L'^{fl), which implies that {un)^^^^ 



converges to in H'i{U), and therefore in Cj,{il.). In particular, u/2 < un < 2 



u 



on VL for N large enough, so that we can assume in our analysis, without loss of 
generality, that F satisfies (6) with q = Q and (8) with r = 2 and s = 0. We also 
deduce from (42) that ujq converges to u in H'^'^{Vt). 

Besides, the unique solution to (23) solves the elliptic equation 
-AV-^ + {V + f{u^) + 2f'iu'y - X) 

f{u^)u^ijj^) u + w - {w,u)l2u, (44) 



n 
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from which wc mfcr that ^u^-u S H'^{Vl) and UTAuN-uHifa < C\\un — u\\l2. Hence, 

1 C 

\\%l)u„-u-^N-4^UN-u\\m < -^||V^«jv-u||h2 < —\\UN - u\\l2. 

Wc therefore deduce from Theorem 1 that 

C 

\\un-u\\h^ < ^„+2^s for s = and s = 1 (45) 
\>^N-\\ < (46) 
^\\uN - u\\jji < E{un) - E(u) < C||u7v - "IIhi- 

From (45) and the inverse inequahty 

which holds true for all s < and all iV > 1, we then obtain using classical arguments 
that 

\\un-u\\h^, < forall0<s<a + 2. (47) 

The estimate (46) is slightly deceptive since, in the case of a linear eigenvalue 
problem (i.e. for ~Au + Vu ~ Xu) the convergence of the eigenvalues goes twice as 
fast as the convergence of the eigenvector in the iJ^-norm. We are going to prove 
that this is also the case for the nonlinear eigenvalue problem under study in this 
section, at least under the assumption that F e cM+^'''"M+'((0, +cx)),]R). 

Let us first come back to (32), which we rewrite as, 

Ajv - A = {{Au - A)(ujv - u), (it AT - u))x'.x + / Wu.«„("Ar - It) (48) 

with 

^ f«)-n^') -,2. , fiu%)~fiu^) 

Wu,UN — ^ Ui^[Un + U) 5 ^ . 

UN — U Upj — 

As ii/2 < u^v < 2w on for N large enough, as u^r converges, hence is uniformly 
bounded, in H^+^{n) and as / e CM+i''^-W+^([||w||ioo/4, 4||m|||^], R), we obtain 
that Wu,UN is uniformly bounded in H^{il) (at least for N large enough). We 
therefore infer from (48) that for N large enough 

|Aiv - A| < C {\\uN - ullffi + \\uN- uWh-") ■ (49) 

Let us now compute the i?~''-norm of the error for < r < cr. Let w E H^{Q). 
Proceeding as in Section 2. wc obtain 

w[uN -u) = {{E"{u) - X){uN - u), n~ .i>w)x',x 



IX'.X 



+ ||uAr-u|||2 / f'{u^)u^1p.u] -]^\\uN - u\\\2 / uw, (50) 

where Hj^ denotes the orthogonal projector on H for the inner 
product. Wc then get from (44) that -0^, is in H^'^{n) and satisfies 

< C\\w\\Hr, (51) 
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for some constant C independent of w. 

Combining (18), (38), (39), (47), (48), (50) and (51), we obtain that there exists a 
constant C G IR+ such that for aU iV € N and all w G H^{Q), 

wiuN-u) < C" {\\uN-u\\l,+N~<^'-+^^\\uN-u\\H^)\\w\\Hr 

(J 

Therefore 

w(un — u) 



C 

\\uN-u\\H-r= sup < (^2) 

for some constant C G IR+ independent of N. Using (47) and (49), we end up with 

C 



\\n-\\ < 



We can summarize the results obtained in this section in the following theorem. 



Theorem 2 Assume that V G H^(Q) for some a > d/2 and that the function F 
satisfies (5)-(7) and is in C['^l+"'"''^^['^l+'^((0, +oo), R). Then {uj^)N^fi converges to 
u in (O) and there exists C G M_|_ such that for all N & N, 

\\un-u\\h.^ < j^^:^ for all -a <s<a + 2 (53) 
IAat-AI < ^ 



'^\\UN - u\\]J^ < E{un) - E{u) < C\\uN-u\\l^. (54) 
//, in addition, F G CW+2'^-W+^((0, +oo), M), then 

\Xn-X\<j^^. (55) 



In order to evaluate the quality of the error bounds obtained in Theorem 2, we have 
performed numerical tests with = (0, 2tt), V{x) = sin(|a;— 7r|/2) and F{t'^) = t'^/2. 
The Fourier coefficients of the potential V are given by 



1 



Vu = -^j—-j. (56) 



4 



from which we deduce that V G H^{Q, 2n) for all a < 3/2. It can be see on Figure 1 
that ||7ijv — \\un — u\\i^2, \\un — u\\fj-i, and \Xn — A| decay respectively as 

jy-2.67 ^ ^-3.67^ ^-4.67 g^j-^^^ j\j--5 ^^]^g reference values for u and A are those obtained 

for N = 65). These results are in good agreement with the upper bounds (53) (for 
s = 1 and s = 0), (52) (for r = 1) and (55), which respectively decay as TV^^-^^'^, 
iV-''-5+^ and for e > arbitrarily small. 
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Figure 1: Numerical errors ||itAr — (+), ||itAr — u||2;,2 (x), ||uAr — (*), and 

\^N — M functions of 2N + 1 (the dimension of Xn) in log scales. 

4 Finite element discretization 

In this section, we consider the problem 

mil^E{v), v€X, J^v^ = lY (57) 

where 

51 is a rectangular brick of R"*, with = 1, 2 or 3, 

^^''^ " ^ X '^"'i' + ^ X + ^ X 

We assume that V E i^(il) and that the function F satisfies (5)- (7), as well as (8) 
for some 1 < r < 2 and < r + s < 3. Throughout this section, we denote by u the 
unique positive solution of (57) and by A the corresponding Lagrange multiplier. 

In the non periodic case considered here, a classical variational approximation of (1) 
is provided by the finite element method. We consider a family of regular triangu- 
lations (7/j)/j of O. This means, in the case when d = 3 for instance, that for each 
h > 0, Tfi is a. collection of tetrahedra such that 

• is the union of all the elements of 7^; 

• the intersection of two different elements of Th is either empty, a vertex, a 
whole edge, or a whole face of both of them; 

• the ratio of the diameter of any element K of to the diameter of its 
inscribed sphere is smaller than a constant independent of h. 

As usual, h denotes the maximum of the diameters hj^, K E Th- The parameter of 
the discretization then is S = h > 0. For each K inTh and each nonnegative integer 
k, we denote by Pfc(if ) the space of the restrictions to K of the polynomials with d 
variables and total degree lower or equal to k. 
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The finite element space X^^k constructed from Th and Pfc(A') is the space of all con- 
tinuous functions on 17 vanishing on 9f2 such that their restrictions to any element 
K of % belong to Pfc(A'). Recall that Xh,k C Hl{Vl) as soon as fc > 1. 

We denote by ttJJ and ii\ j, the orthogonal projectors on Xh^k for the and 
inner products respectively. The following estimates are classical (see e.g. [8]): 
there exists C £ M+ such that for all r G N such that 1 < r < A; + 1, 

y<t>eH-in)nH^in), u-7rl,(f>\\„i <ch^-'u\\H''. (58) 

Let Uh,k be a solution to the variational problem 

,,2 



\\uh,i - u\\h^ 


< 


Ch 


\\uh,i - u\\l2 


< 


Ch^ 




< 


Ch^ 


u\\hi < E{uh,i) - E(u) 


< 


Ch^ 



inf <yE{vh,k), Vh.k e Xh,k, J U/i^fe = 1 

such that {uh,k,u)L2 > 0. In this setting, we obtain the following a priori error 
estimates. 

Theorem 3 Assume that V G L^(il) and that the function F satisfies (5), (6) for 
q = 1, (7), and (8) for some 1 < r < 2 and < 7' + s < 3. Then there exists /iq > 
and C G R+ such that for all < h < Hq, 

(59) 
(60) 
(61) 

l\\uhA-u\\j,, <E{uh.i)-E{u) < Ch^. (62) 

// in addition, V G H^(fl), F satisfies (8) for r = 2 and is such that F G 
C3((0,+oo),R) and F"{t)t^/'^ and F"'(t)t^/'^ are bounded in the vicinity ofO, then 
there exists ho > and C G M+ such that for all < h < ho, 

\\uh,2-u\\Hi < Ch-" (63) 
\\Uh,2~u\\L2 < Ch^ (64) 
|A„,2-A| < Ch^ (65) 

^\\uh,2-u\\jj, <E{uk.2)-E{u) < Ch\ (66) 

Proof As is a rectangular brick, V satisfies (4) and F satisfies (5)-(7), we have 
u G H^{^1). We then use the fact that ipu^ ^-u is solution to 

-Ai^u,^,.^ + iV + f{u') + 2f'{u^)u' ~ \)4>u,,,-u 

= 2 f'{u^)u^%l;uu,k-u^ + ("/j^fe - w) - {uh.k - u, u)l2U, 

to establish that ipuj^ t-M G H^{fl) n H^{Q) and that 

< C\\uh,k - u\\l^- (67) 

for some constant C independent of h and fc. The estimates (59)-(62) then are 
directly consequences of Theorem 1, (30), (58) and (67). 

Under the additional assumptions that V G H^{n), we obtain by standard elliptic 
regularity arguments that u G iJ'^(ri). The and estimates (63) and (64) 
immediately follows from Theorem 1, (30), (58) and (67). We also have 

\X2,h - A| < C/i^ 
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for a constant C independent of h. In order to prove (65), we proceed as in Section 3. 
We start from the equality 



A2,/i - A = {{Au - X){u2,h - {U2,h - u))x',X + W (li2,/i - u) 

Jn 

where 

U2,h - U 

We now claim that Uh,2 converges to u in L°°{n) when h goes to zero. To establish 
this result, we first remark that 

\\uh,2 - uIIl- < \\uh,2 ~^h,2u\\L-' + \\^h,2U - u\\l'^, 

where Th,2 is the interpolation projector on Xh^2- As u £ H^{il) ^ C^(il), we have 

lim \\2h 2U — — 0. 

On the other hand, using the inverse inequality 

3C e IR+ s.t. VO < /i < ha, e X,,,2, hh^L^ < Cp{h)\\vh,2\\H^ , 

with p{h) = 1 if = 1, p{h) = 1+lnhif d = 2 and p{h) = h-^/^ if d = 3 (see [8] 
for instance), we obtain 

\\uh,2 -Ih,2u\\L'== < Cp{h)\\uh,2 -Ih.2u\\Hi 

< Cp{h) {\\uh,2 - u\\h^ + ||u - Jh,2u|lHi) 

< C'p(h)h^ — > 0. 

/wO+ 

Hence the announced result. This implies in particular that w'* is bounded in H^{il), 
uniformly in h. Consequently, there exists C E M+ such that for all < h < Hq, 

|Ah,2 - A| < C (||uh,2 - + \\uh,2 - "llff-i) • (68) 
To estimate the iJ^^-norm of Uh 2 — u, we write that for all w € 7?g(f2), 

w{uh^2-u) = {{E"{u) - X){uh^2 - u),'k]^^^^^^±iPw)x',x 

+ {{E"{u) - X){uh,2 - u), Vu- - T^x^^^nu^'^-'-w)x'.x 

+ \\Uk,2-u\\l2 I /'(m^)u^V'^ - ^||w/i,2 -U|li2 / UW, 

Ja ^ Jn 

where ipy^ is solution to 

-AV^„, + {V + f{u^) + 2f'{u^)u^ - A)V'„, 

f'{u'^)u^'ipw) u + w - {w,u)l2u, (69) 
n J 

and where tt^^ jHu^ denotes the orthogonal projector on Xh_2 H for the 
inner product. Using the assumptions that V € H^{i}), F g C'^((0, +oo), E), and 
F"{t)t^/^ and F"'{t)t^/'^ are bounded in the vicinity of 0, we deduce from (69) that 
i/'u, is in H^iyi) and that there exists C £ R+ such that for all w £ Ho{^) and all 
<h<ho, 

W^^Wh^ <C\\w\\m. 
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We therefore obtain the inequaUty 

U^~7tI^^J\hi < Ch^\\w\\m, (70) 
where the constant C is independent of h. 

Putting together (8) (for r = 2), (18), (38), (39), (58), (63), (64) and (70), we get 

IK. --Ik— sup k^f!;^<ch\ 

Combining with (63) and (68), we end up with (65). Lastly, we deduce (66) from 
the equahty 

E{uh.2)-E{u) = ]^{{Au-X)[Uh,2-u),{uh,2~u))x',X 



Taylor expanding the integrand and exploiting the boundedness of the function 
F"{t)t^/'^ in the vicinity of 0. □ 

Numerical results for the case when Q, = (0, tt)^, V{xi, X2) = x\~\- and F{t^) ~ 
t^/2 are reported on Figure 2. The agreement with the error estimates obtained in 
Theorem 3 is good for the Pi approximation and excellent for the P2 approximation. 



5 The effect of numerical integration 

Let us now address one further consideration that is related to the practical imple- 
mentation of the method, and more precisely to the numerical integration of the 
nonlinear term. For simplicity, we focus on the case when A~\. 

From a practical viewpoint, the solution (u^-, A5) to the nonlinear eigenvalue prob- 
lem (15) can be computed iteratively, using for instance the optimal damping algo- 
rithm [4, 2, 7]. At the p'^ iteration {p > 1), the ground state (w^, A^) e x M of 
some linear, finite dimensional, eigenvalue problem of the form 

(v^-Vvs+j (V + fip'f')) ^vs = \l [ ^vs, yvs e Xs, (71) 
Jn Jn ^ ^ Jo 

has to be computed. In the optimal damping algorithm, the density is a 

convex linear combination of the densites pg = |w'P, for < < p — 1. Solving 
(71) amounts to finding the lowest eigenelement of the matrix with entries 



iJP / V0fe-V0z+ V<t>k<t>i+ f{pri<t>k<t>u (72) 
Jn Jn Jn 

where ('/'fc)i<j,<(jini(X5) stands for the canonical basis of Xs- 

In order to evaluate the last two terms of the right-hand side of (72), numerical 
integration has to be resorted to. In the finite element approximation of (57), it 
is generally made use of a numerical quadrature formula over each triangle (2D) 
or tetrahedron (3D) based on Gauss points. In the Fourier approximation of the 
periodic problem (40), the terms 



/ Vekei and / /(p^ ^)ek< 
Jn Jn 
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10' 




Figure 2: Errors \\uh.k - (+), \\uh.k ~ u\\l2 (x) and \Xh,k ~ A| (*) for the Pi 

{k = 1, top) and P2 (fc = 2, bottom) approximations as a function of h in log scales. 



which are in fact, up to a multiplicative constant, the (fc — Z)"^ Fourier coefficients of 
V and f{pg^^) respectively, are evaluated by Fast Fourier Transform (FFT), using 
an integration grid which may be different from the natural discretization grid 

{ (2^^'^' 2^^"^' )' ° '"^1 

associated with Xjy. This raises the question of the influence of the numerical 
integration on the convergence results obtained in Theorems 1, 2 and 3. 

Remark 4 In the case of the periodic problem considered in Section 3 and when 
F{t) = ct^ for some c > 0, the last term of the right-hand side of (72) can he 
computed exactly (up to round-off errors) by means of a Fast Fourier Transform 
(FFT) on an integration grid twice as fine as the discretization grid. This is due to 
the fact that the function ^~^'ekei belongs to the space Span{en \ < 4N}. An 
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analogous property is used in the evaluation of the Coulomb term in the numerical 
simulation of the Kohn-Sham equations for periodic systems. 



In the sequel, we focus on the simple case when d = 1, = (0, 2tt), X = H^{0, 27r), 
and 

1 



E{v) 











Jo 




4 70 



27r 

L,|4 



with V e H^{Q,2tt) for some <t > 1/2. More difficult cases will be addressed 
elsewhere [5]. 

In view of Remark 4, we consider an integration grid 

27r „ , , r 27r 47r 27r(7V„ - 1) 
— Zn[0,27r)= 0,— 



with Ng > AN + 1 for which we have 

r2n 



^ ^6^zn[o,2^) 

and for all p G X2N^ 

V|fcUI|<iV, rpej:ei^^ ^ PWe"'^''')'- = (73) 

-zn[o,2.) 

where is the (fc — Z)'^ coefficient of the discrete Fourier transform of p. Recall 
that if = X]gGZ'?^9^9 ^ C'^(0,27r), the discrete Fourier transform of 4> is the 
A'^gZ-periodic sequence defined by 

^ '■6-^Zn[0,27r) 

We now introduce the subspaces Wm for M E N* such that Wm = -'^(j\f-i)/2 if M is 
odd and M^m = ^m/2-i © C(eM/2 + e-M/2) is is even (note that dim(Wj\/) = M 
for all M G N*). It is then possible to define an interpolation projector Xjv from 
C^(0,27r) onto VKjVg by 

Va; G I^Zn [0,27r), [liv,(0)](x) = (b{x). 

The expansion of (0) in the canonical basis of is given by 

(2^)1/2 ^ (TV, odd), 

\g\<{N, -l)/2 

(2.)i/2 ^ ^ + (2.) V^0FPT 

|g|<JVg/2-l ^ ^ 



{Ng even). 



Under the condition that Ng > AN + 1, the following property holds: for all e 
C0(0,27r), 



27r 



v|fc|,|;|<iv, / ImM^^^ = ^VJ 
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l'2-iT 


p2TT 






/ f^N 




Jo 


Jo 



It is therefore possible, in the particular case considered here, to efSciently evaluate 
the entries of the matrix using the formula 

RP - 
^ki ■- 

- (74) 

and resorting to Fast Fourier Transform (FFT) algorithms to compute the discrete 
Fourier transforms. Note that only the second term is computed approximatively. 
The third term is computed exactly since, at each iteration, belongs to X2N 
(see Eq. (73)). Of course, this situation is specific to the nonlinearity F{t) = 
considered here. 

Using the approximation formula (74) amounts to replace the original problem 

inf \e{vn). vn e Xn. \vNf = l| , (75) 
with the approximate problem 

inf Ien^vn), vn e Xn, f \vn\^ = l] , (76) 



where 

Let us denote by un a solution of (75) such that (ujv,u)i2 > and by UN,Ng 
a solution to (76) such that (mtv,^,, ■")^2 > 0. It is easy to check that UM,Ng is 
bounded in H^{0, 2tt) uniformly in and Ng. 

Besides, we know from Theorem 2 that (ujv)7VeN converges to u in H^(0, 27r), hence 
in (2, tt) , when goes to infinity. This implies that the sequence (yl„ — ) at^n 
converges to in operator norm. Consequently, for all TV large enough and all Ng 
such that Ng>iN+l, 

'^\\uN,Ng - UN\\jji < E{un,nJ - E{un) 

< {un,nJ - Epfg (un) 

'''\v-lNAV)){\uN,Nf-\uN\') 

Jo 

< C\\I12n{V ~ lNAV))\\L4uN,Ng ~ UN\\m, 

where we have used the fact that {\uN^]y^\'^ — |u7vP) S ^2JV- Therefore, 

\\uN,Ng-UN\\H^ <C\\Il2N{V-lNg{VmL^. (77) 

for a constant C independent of and Ng. Likewise, 

AjV,JVg — = ((^ujv ^ ^N){uN,Ng ~ Un), (wjv.ATg — Un))x',X 

+ f (V~lN{V))\uN,Nf 

Jo 

l'2-K 
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from which wc deduce, using (77), 



\Xn,n,~Xn\ <C\\U2N{V-lNJVm 



L2- 



An error analysis of the interpolation operator XjVg is given in [6]: for all non- 
negative real numbers < r < s with s > 1/2 (for d ~ I), 

\\^-lN,{^)\\Hr < -^II^IU^, V^ei/^(0,2^). 

Thus, 

\\n2N{V-lNAV))\\L^ < \\V~InAV)\\l^ < (78) 
and the above inequality provides the following estimates: 

\\un,n,-u\\hi < C{N-^-'+N-^) (79) 
\\um,n,-u\\l2 < C {N-'^-^ + (80) 
\Xn,n,-X\ < C (iV-2-"2 ^ iv--) , (81) 

for a constant C independent of N and Ng. The first component of the error bound 
(79) corresponds to the error \\un — u\\h'^ while the second component corresponds 
to the numerical integration error \\un^n —un\\h^ (the same remark applies to the 
error bounds (80) and (81)). 

It is classical that for the norm \\ip — Tj^j^lp\\h^- for r < is in general of the same 
order of magnitude as \\ip—2]^^ip\\]^2. As the existence of better estimates in negative 
norms is a corner stone in the derivation of the improvement of the error estimate 
(46) for the eigenvalues (doubling of the convergence rate), we expect that the 
eigenvalue approximation will be dramatically polluted by the use of the numerical 
integration formula. 

This can be checked numerically. Considering again the one-dimensional example 
used in Section 3 (f^ = (0,27r), V{x) =^ sin(|x - 7r|/2), F{t) = we have 

computed for 4 < iV < 30 and Ng = 2^ with 7 < p < 15, the errors Wun^n^^ — 
\\uN,Ng — llwAT^ATg — and lAAr^Ar^ — A|. On Figure 3, these quantities are 

plotted as functions of 2N + 1 (the dimension of X^), for various values of Ng. 

The non-monotonicity of the curve TV i— > IAat^at^ — A| originates from the fact that 
Aw.Wg — A can be positive or negative depending on the values of N and Ng. 



The numerical errors Wun.n^-uWh^, llww.Wg-MlU^, \\un^n^-u\\h-i , and |AAr,Arg-A|, 
for N = 30, as functions of Ng (in log scales) are plotted on Figure 4. When Ng goes 
to infinity, the sequences log^, \\uN,Ng - u\\m, log^o \\uN,Ng - logio W^n^n^ - 

uWh-i, and log^o IAat.at^ - A| converge to log^o ||uAr - log^p ||ujv - "lU^, 

logj^Q \\un — and logj^Q \Xn — A| respectively. For smaller values of iV^, the 

numerical integration error dominates and these functions all decay linearly with 
log^o A'g with a slope very close to —2. For fixed N, the upper bounds (79)-(81) 
also decay linearly with logj^g but with a slope equal to —1.5. To obtain sharper 
upper bounds for the numerical integration error, we need to replace (78) with 
a sharper estimate of ||n27v(^ — ^Ng{V))\\L2, which is possible for the particular 
example under consideration here. Indeed, remarking that under the condition 

Ng >47V+1, 

\\U2N{V-lNg{V))h'^ = I 

\\9\<2N 



E 



g+kNg 



1/^ 
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3. 10 4. 10 5. 10 6. 10 





3. 10 4. 10 5. 10 6. 10 



3.10 4.10 S. 10 6.10 7.10 



Figure 3: Numerical errors \\un^n^ — u\\f{i (top left), \\uN,Ng — u\\l^ (top right), 
||itAf,Afg — wjl//-! (bottom left), and |AArjVg — M (bottom right), as functions of 
2N + 1 (the dimension of Xn), for Ng = 128 (red), Ng = 256 (green), Ng = 512 
(cyan), Ng = 1024 (gold), Ng = 2048 (magenta), Ng = 4096 (pink), Ng = 8192 
(black), Ng = 16384 (blue), Ng = 32768 (light blue). 



we can, using (56), show that 



\\Il2NiV ~Jn,{V))\\l2 < 



CN^ 
N^ ' 



for a constant C independent of N and Ng. We deduce that for this specific example 

\\uN,N,~uU2 < c(^N-y^ + N^/^N-'') 



\Xn.n,-X\ < c[n-'^' + N'/'N-') 



Acknowledgements 

This work was done while E.G. was visiting the Division of Applied Mathematics 
of Brown University, whose support is gratefully acknowledged. The authors also 
thank Jean- Yves Chemin and Didier Smets for fruitful discussions, and Claude Le 
Bris for valuable comments on a preliminary version of this work. 



23 



-9 - 
10' 



Figure 4: Numerical errors lluAr^Ar^ - m||_h-i (x), \\un,n^ -uWl-^ (+), \\uN^Ng-u\\ 
{*), and |Ajv,jVg — M (°)' foi' ^ = 30, as functions of Ng (in log scales). 
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6 Appendix: properties of the ground state 

The mathematical properties of the minimization problems (1) and (9) which are 
useful for the numerical analysis reported in this article are gathered in the following 
lemma. 

Recall that = 1. 2 or 3. 



Lemma 2 Under assumptions (2)-(6), (9) has a unique minimizer po and (1) has 
exactly two minimizers u = and —u. The function u is solution to the nonlinear 
eigenvalue problem (11) for some A G M. Besides, u S C°'"(r2) for some < a < 1, 
u > Q in n, and A is the lowest eigenvalue of and is non-degenerate. 



Proof As A is uniformly bounded and coercive on Vi and V € L'^(r2) for some 
q > max(l,c?/2), v i-^ a{v,v) is a quadratic form on X, bounded from below on 
the set {v € X \ \\v\\l2 — 1}. Replacing a[v, v) with a[v, w) + C||w||^2 and F{t) with 
F [t] ~ F [0] — tF' [0) does not change the minimizers of (1) and (9). We can therefore 
assume, without loss of generality, that 

Vw G X, a{v, v) > \\v\\l2 and F(0) = F'{0) = 0. (82) 

It then follows from (6) and (82) that < F{v^) < C{v^ + v^). As X ^ L^{^), 
E{v) is finite for all f G X, / > — oo and the minimizing sequences of (1) are 
bounded in X. Let (w„)ngN be a minimizing sequence of (1). Using the fact that 
X is compactly embedded in L^(f2), we can extract from (i'„)neN a subsequence 
(w„^)fegN which converges weakly in X, strongly in L^(ri) and almost everywhere 
in f2 to some u € X. As llwri^llLa — 1 and E{vn^) [ I, we obtain \\u\\]^2 = 1 and 
E{u) < I {E is convex and strongly continuous, hence weakly l.s.c., on X). Hence 
M is a minimizer of (1). As |m| G X, |||M|||i2 = 1 and i?(|M|) = E{u), we can assume 
without loss of generality that u > 0. Assumptions (2)-(6) imply that E is on 
X and that E'{u) = A^u. It follows that u is solution to (10) for some A G M. By 
elliptic regularity arguments [10], we get u G C°'"(0) for some < a < 1. We also 
have u > in il; this is a consequence of the Harnack inequality [12]. Making the 
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change of variable p = u^, it is easily seen that if w is a minimizer of (1), then is a 
minimizer of (9), and that, conversely, if p is a minimizer of (9), then y/p and 
are minimizers of (1). Besides, the functional E is strictly convex on the convex set 
{p > I 7p e X, p = 1}. Therefore po = V? is the unique minimizer of (9) and 
u and ~u are the only minimizers of (1). 

It is easy to see that A„ is bounded below and has a compact resolvent. It therefore 
possesses a lowest eigenvalue Aq, which, according to the min-max principle, satisfies 

Ao = inf ( / (AVw) • + / (y ^ v(^X, f = l\ . (83) 

Let vq be a normalized eigenvector of associated with Ao- Clearly, vq is a 
minimizer of (83) and so is |uo|. Therefore, |wo| is solution to the Euler equation 
^ul^ol — '^ol'^^ol- Using again elliptic regularity arguments and the Harnack in- 
equality, we obtain that |uo| € C°'"(il) for some < a < 1 and that \va\ > on 17. 
This implies that either va ~ \va\ > in 17 or wo = ~K'o| < in il. In particular 
(m, vq)i2 7^ 0. Consequently, A = Aq and A is a simple eigenvalue of A^. □ 



Let us finally prove that A is also the ground state eigenvalue of the nonlinear 
eigenvalue problem 

search (p, w) e K x X such that 

AyV = pw (84) 

ll«IU= = i, 

in the following sense: if (p, v) is solution to (84) then either p > A or p = A and 
V = ±u. 

To see this, let us consider a solution (p, v) € MxX to (84) and denote by ii; = |?;|— it. 
As for It, we infer from elliptic regularity arguments [10] that v E C°'"(ri). We have 
= ||u||i2 = 1. Therefore, if w < in f2, then = u, which yields v = ±u and 
p = A. Otherwise, there exists xq E such that w{xo) > 0, and, up to replacing v 
with —V, we can consider that the function w = v — u is such that w{xo) > 0. The 
function w is in X D C'^'"(n) and satisfies 

[Au - X)w + /(^') -/("') + = _ x)v. (85) 
V — u 

Let u) = {x G O I w{x) > 0} = {x Eil \ v{x) > u{x)} and w+ = max(w,0). As 
w+ e X, we deduce from (85) that 

({Au - \)w+,w+)x'.x + / v{u + v)w'^ = - X) / vw. 

Jul ^ ^ Juj 

The left hand side of the above equality is positive and vw > 0. Therefore, p > A. 
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